Photoassociation and coherent transient dynamics in the interaction of ultracold 
rubidium atoms with shaped femtosecond pulses - II. Theory 
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Photoassociation of ultracold rubidium atoms with femtosecond laser pulses is studied theoret- 
ically. The spectrum of the pulses is cut off in order to suppress pulse amplitude at and close to 
the atomic resonance frequency. This leads to long tails of the laser pulse as a function of time 
giving rise to coherent transients in the photoassociation dynamics. They are studied as a function 
of cutoff position and chirp of the pulse. Molecule formation in the electronically excited state is 
attributed to off-resonant excitation in the strong-field regime. 
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I. INTRODUCTION 

Cooling, trapping and manipulation of atoms and 
molecules in the ultracold regime (T < 100 //K) repre- 
sents one of the most active research fields in contempo- 
rary atomic, molecular, and optical physics. Photoas- 
sociation is the formation of molecules when two col- 
liding atoms are excited by laser light into bound vi- 
brational levels of an electronically excited state. Us- 
ing continuous- wave lasers, photoassociation of ultracold 
atoms can be employed to study collisional properties of 
ultracold atoms, in particular long-range potential en- 
ergy curves [l|. If followed by spontaneous emission, it 
also serves to form ultracold molecules in their electronic 
ground state [1]. 

The use of short laser pulses for photoassociation of 
ultracold atoms has been suggested in theory work for 
about a decade I, i, i, 1 0, ill, [13 . Initially, the mam 
goal was to find a photoassociation scheme that would al- 
low for larger formation rates of excited state molecules 
[3, [3| • In order to coherently produce ultracold molecules 
in their electronic ground state, the pump-dump scenario 
was adapted to photoassociation j6|,[3, [3: A pump pulse 
creates a vibrational wavepacket on the electronically ex- 
cited state which travels toward its inner turning point 
where it is sent to the ground state by the dump pulse. 
Other coherent control schemes such as adiabatic pas- 
sage [§] and pulse optimization using genetic algorithms 
were investigated as well [To|. Coherent photoassocia- 
tion using short laser pulses is closely linked to the ulti- 
mate goal of ultracold molecule formation, the stabiliza- 
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tion of the molecules into their vibrational ground state 
[ill [1^ [13I [l3 | . This has been discussed in theory work 
in the context of optimal control [is'l and coherent pulse 
accumulation [l6j . 

In the theoretical studies on short-pulse photoassocia- 
tion, pulses with a spectral bandwidth of a few wavenum- 
bers were suggested. This choice is motivated by the 
requirement of exciting a narrow band of transitions 
with large free-bound Franck-Condon factors close to 
the atomic resonance without exciting atomic transi- 
tions. This spectral bandwidth corresponds to transform- 
limited pulse durations of a few picoseconds. However, 
pulse shaping capabilities for such pulses have yet to be 
developed. 

Femtosecond pulses on the other hand can be shaped 
with a number of techniques, but they have a very broad 
spectrum, addressing both molecular and atomic transi- 
tions. Since the probability of atomic transitions is sev- 
eral orders of magnitude larger than that of photoasso- 
ciation transitions, the pulse spectral amplitude at the 
atomic resonance frequency needs to be completely sup- 
pressed in order to record any molecular signal. This 
is most easily achieved by putting a knife edge into the 
Fourier plane of the pulse shaper [l3| ■ The correspond- 
ing sharp cutoff in the spectrum of the femtosecond pulse 
leads to long tails of the pulse as a function of time and 
induces coherent transient dynamics [l^, [3| ■ The tran- 
sients are reflected in oscillations of the molecular pop- 
ulation and can be measured by ionization detection of 
the excited state molecules [l^l • Photoassociation of ul- 
tracold atoms with shaped femtosecond laser pulses has 
thus been demonstrated recently in a break-through ex- 
periment [20| . 

This paper presents a detailed explanation of the ex- 
perimental results obtained by femtosecond photoasso- 
ciation of ultracold rubidium atoms with shaped pulses 
based on theoretical calculations. The details of the ex- 



perimental setup as well as the ma jor part of the ex- 
perimental results is found in Ref. [2l|. Figure [1] de- 
picts the proposed photoassociation mechanism includ- 
ing the potential energy curves considered in the calcu- 
lations. A first laser field (the pump pulse, indicated by 
the red arrow) excites two atoms from their ground state 
(5s-|-5s) into long-range bound molecular levels of an elec- 
tronically excited state (correlated to the 5s-|-5pi/2 or 
5s-|-5p3/2 asymptotes). The initial scattering state with a 
relative scattering energy corresponding to the trap tem- 
perature of 100 fiK is shown on the ground state. The 
excited state wavefunction after the pulse is over as well 
as one exemplary excited state vibrational wavefunction 
are displayed in the first excited state of Fig. [T] In order 
to detect the excited state molecules, a second laser field 
(the probe pulse, indicated by the green arrow) transfers 
them to the ionic state. 

The paper is organized as follows: Section [IT] briefly re- 
views the theoretical description of two rubidium atoms 
excited by a photoassociation laser pulse. The photoas- 
sociation mechanism based on strong-field off-resonant 
excitation is described in detail in Sec. IIIII while the co- 
herent transient dynamics are studied in Sec. II VI Theo- 
retical and experimental results are compared in Sec. IVl 
and Sec. IVII concludes. 

II. THEORETICAL APPROACH 

In order to gain detailed insight into the mechanism 
of femtosecond photoassociation, the pump excitation 
step was simulated in quantum dynamical calculations 
of the light-molecule interaction. The pump pulse trans- 
fers population between the ground and first excited elec- 
tronic states. The probe step was not explicitly taken 
into account. Rather, it was assumed that any popula- 
tion of the first excited state is ionized alike. 

In a two-channel picture, the Hamiltonian describing 
two rubidium atoms subject to excitation by a laser field 
in the dipole and rotating-wave approximations is given 

by 

\-fiE{tr T + Ve{R)-A)- ^ ^ 

Here, R denotes the internuclear distance, T the kinetic 
energy operator and Vg/e the ground and excited state 
potential energy curves. The parameter A=h{ujo — ujl) 
represents the detuning of the pump pulse central fre- 
quency relative to the 5s -I- and 5s -I- 5p3/2 asymp- 
totes, respectively. The off-diagonal elements of the 
Hamiltonian describe the electric dipole coupling with /t 
denoting the transition dipole and E{t) the electric field 
of the pump pulse. The i?-dcpcndcncc of the transition 
dipole was neglected which is justified for excitation at 
long range. 

A two-channel picture assumes that no couplings be- 
tween electronic states, for example due to spin-orbit in- 
teraction, play a role on the timescale of the dynamics. 
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FIG. 1: (color online) Scheme of the femtosecond pump-probe 
sequence involving the Rb2 potential energy curves of the 
electronic ground state, the attractive potentials of the first 
electronically excited state correlating to the 5s-|-5pi/2 and 
5s-|~5p3/2 asymptotes, and the cationic state of RbJ. The 
initial scattering state, the wavefunction amplitude after the 
pump pulse and one sample excited state vibrational wave- 
function (v = 14 of the 0^(5s + 5p3/2) state) are also depicted. 



It is then possible to consider pairwise combinations of 
ground and excited electronic states in Eq. ([T]). An analy- 
sis of the dynamics has confirmed that this assumption is 
justified: The dynamics are purely long-range and dom- 
inated by electronic transitions due to the laser field. 

The potential energy curves are taken from Ref. (2^ 
and matched to the long-range dispersion potentials, 
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{C3/R^+)Ce/R^ + Cs/R^ + Cw/R^^ with the C„ coeffi- 
cients found in Refs. [H, [13]. Potential energy curves in 
the Hund's case (c) representation were obtained by ap- 
proximating the spin-orbit couplings by their asymptotic 
values and subsequent diagonalization, see e.g. Ref. ^2^. 
Of all potential curves correlating to the 5s -I- 5pi/2 
and 5s -I- 5p3/2 asymptotes, those with attractive char- 
acter and dipole-allowed transitions were taken into ac- 
count (Og", Ig and 0+ correlating to the 5s-|-5pi/2 asymp- 
tote and Oj, Ig, 0+ and correlating to the 5s-|-5p3/2 
asymptote). Moreover simulations were also performed 
for two repulsive potential energy curves, those of the 
lu(5s -I- 5^3/2) and 2g(5s -I- 5^*3/2) states. Note that the 
excited state potentials differ by their effective C3 coeffi- 
cients and strengths of transition dipole moment. 

The Hamiltonian, Eg. ([1]), is represented on a mapped 
Fourier grid [H, [23|. Vibrational wavefunctions such 
as the one for w = 14 of the 0~(5s -I- 5^3/2) excited 
state shown in Fig. [1] are obtained by diagonalizing the 
Hamiltonian with E{t) set to zero. The photoassocia- 
tion dynamics were studied by numerically solving the 
time-dependent Schrodinger equation (TDSE), 

^f^^(-^,{R;t)\_^(^g{R;t)\ 

using the Chebyshev propagator [l!]. ^g/^{R]t) de- 
notes the ground (excited) state component of the time- 
dependent wavefunction. The initial state is taken to 
be the scattering state with scattering energy closest to 
Escatt/ks = Ttrap whcrc ks denotes the Boltzmann con- 
stant and Ttrap the temperature of the trap. In principle 
an average over all thermally populated scattering states 
needs to be performed @. However, since the dynam- 
ics are dominated by electronic transitions, the different 
nodal structures of the scattering states at very large 
distances do not play a role and simulations for a single 
initial state are sufficient. 

The pulses entering Eqs. ^ and ([2|) were obtained sim- 
ilar to the experimental procedure. In the experiment an 
optical low pass filter was employed in the Fourier plane 
of the zero dispersion compressor to manipulate the spec- 
tral shape of the pump pulse close to the Dl or D2 atomic 
line in such a way, that all frequencies above the chosen 
atomic resonance, the respective resonance frequency it- 
self, and some frequencies below this particular line were 
cut off. This is done in order to prevent the loss of the 
captured atoms due to resonant three-photon ionization 
and photon pressure [l3|. In the case of excitation near 
the D2 line, perturbing frequencies near the Dl atomic 
resonance due to the large spectral bandwith of the fem- 
tosecond pulse were also removed from the pump pulse 
spectrum. In order to model those experimental pump 
pulses, a Gaussian transform-limited laser field is sub- 
jected to one or two cuts of spectral amplitude in the 
frequency domain. The sharp cuts lead to very long tran- 
sients of the pulse in the time-domain representation, cf. 
Fig. m i.e. the temporal behavior strongly deviates from 
a transform limited profile. The pulse has long tails. 



reaching out to a few picoseconds before and after the 
maximum as shown in Fig. [51 If two cuts are applied, 
the tails show a beat pattern that depends on the spe- 
cific cut positions, cf. Fig. [2Kb). Within the tails, the 
electric field oscillates with a frequency that matches the 
spectral cutoff frequency, cf. Fig. [3) The influence of the 
specific shape of the pump pulses in time and frequency 
domain on the dynamics will be discussed in Sec. lIVI The 
parameters defining the pulse are the central wavelength, 
spectral width, energy, linear chirp rate and spectral cut- 
off positions. 

The time-dependent population of the first electroni- 
cally excited state is calculated as weighted sum of the 
possible transitions from the singlet ground and lowest 
triplet state with the ratio 1:3. In order to take the fi- 
nite duration of the probe pulse into account, the mea- 
sured ion signal was obtained by convoluting the time- 
dependent excited state population with a Gaussian of 
500 fs full- width at half maximum (FWHM)[3l]. 

III. PHOTOASSOCIATION 

The photoassociation dynamics is illustrated in Fig. [1] 
by the nuclear wavefunctions, ^'g(i?;t) and ^'e(i?;t), 
in the ground and excited state potentials, when the 
pulse is over, at about t = 18 ps after the pump pulse 
maximum. The aS+(5s -I- 5s) lowest triplet state and 
Og (5s -I- 57)3/2) excited state were employed in this exam- 
ple. Pump pulse parameters like pulse energy and cutoff 
position correspond to those used in the actual experi- 
ment. Specifically, the cutoff positions for the example 
displayed in Fig. [T] and Fig. [3Kb), were — IScm"-'^ and 
— 218.3 cm~^ below the D2 resonance, the central wave- 
length Al = 785 nm, spectral bandwidth AAl = 25 nm, 
laser power P = 30 mW, repetition rate 100 kHz and 
beam radius = 300 /um, corresponding to a peak field 
amplitude of 99MV/m. Inspection of Fig. [1] reveals that 
no nuclear dynamics has taken place on the time scale of 
the pulse including its long tails: The excited state wave 
packet, created by the pump pulse, completely reflects 
the nodal structure of the initial scattering state at long 
range. This is not surprising since the vibrational peri- 
ods of weakly bound excited state levels are on the order 
of several tens of picoseconds, while the ground state dy- 
namics evolves even more slowly. 

The excited state wave packet at a certain time can be 
projected onto the vibrational eigenfunctions of the cor- 
responding excited state potential. This yields the vibra- 
tional distribution and proves that photoassociation, i.e. 
the excitation of free colliding atoms into bound levels of 
an excited state, has taken place. The vibrational dis- 
tribution after the pulse is shown for the 0+(5s + 5pi/2) 
state in Fig. [SKa) and for the 0~(5s -|- 5^3/2) state in 
Fig. [SKb). The cutoff position was chosen to be 15cm~^ 
below the respective atomic resonance, i.e. below the 
Dl line frequency in Fig. ^a,) and below the D2 line 
frequency in Fig. [3Kb). In Fig. [3Ka), the vibrational dis- 
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FIG. 2: (color online) Temporal field envelope of cut pump pulse spectrum (shown in the insets) below the Dl (a) and D2 (b) 
atomic resonance. Tails to both sides of the peak are caused by the spectral cutoff. 
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FIG. 3: (color online) Calculated population distribution in the excited 0^ (5s + 5pi/2) (a) and Og (5s + 5^3/2) (b) states. The 
solid line displays the spectral pump pulse amplitude where the transition energies from the initial ground state scattering state 
to the bound excited state vibrational levels were matched to the pump pulse spectrum. 



tribution is characterized by a maximum near v' = 400. 
These vibrational levels have large bond lengths on the 
order of 1000 oq and low binding energies of about 1 GHz 
or 0.03 cm^^. They lie energetically above the spectral 
cutoff of the pulse. Their population therefore necessar- 
ily arises from off-resonant excitation. This is caused by 
the high peak intensity of the pump pulse. Only a small 
fraction of the total excited state population resides in vi- 
brational levels that are resonant within the pump pulse 
spectrum. For a spectral cutoff detuned by -15 cm~^ from 
the atomic line, resonant excitation occurs for vibrational 
levels below v' = 290, cf. inset of Fig. [Sja). Compared 
to the weakly bound, long range levels around v' = 400, 
their population is suppressed by two orders of magni- 
tude. This is due to much smaller free-bound Franck- 
Condon factors. The distribution among these levels re- 
produces the weak-field result, i.e. it is proportional to 
the Franck-Condon factors multiplied by the spectral en- 



velope of the pump pulse [13, HI]. The apparant con- 
tradiction of producing within one calculation weak-field 
(below v' = 290) and strong-field (near v' = 400) results 
is easily resolved by noting that the decisive quantity is 
the Rabi frequency, i.e. the product of transition dipole 
and electric field. Since the Franck-Condon factors vary 
over several orders of magnitude, qualitatively different 
results are obtained for different regions of the excited 
state vibrational spectrum. 

A pattern similar to that of the Dl line results is ob- 
served for the D2 line, cf. Fig. [3l^a) and (b). In this 
case, the excited state potential, the purely long range 
0~(5s -|- 5^3/2) state, is only about 28cm~^ deep and 
resonant excitation occurs for vibrational levels below 
v' = 15. Due to the shallowness of the potential, only 
a small part of the pulse spectrum serves for excitation, 
cf. the inset of Fig. [2l^b). If an excited state potential 
with a deep minimum such as that of the lg(5s -I- 5p3/2 
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state is employed (data not shown), a picture equiva- 
lent to Fig. [3Ja) emerges. In both cases the vibrational 
spectrum is dominated by levels which are off-resonantly 
excited, above v' = 90 for the 0~(5s + 6^3/2) state (with 
binding energies of less than 0.05 cm~^ or 1.5 GHz). It 
is instructive to compare the amplitudes of the excited 
state vibrational level with v' = 14 and the excited state 
wave packet after the pulse, shown in Fig. [U to the vi- 
brational distributions of Fig. [3)Jb) . The short range part 
of the wave packet at internuclear distances R < 50 oq 
corresponds to vibrational levels with v' < 14. The to- 
tal excited state population is dominated by the weakly 
bound levels and electronically excited continuum states. 
Population of the latter, however, corresponds to excita- 
tion of the initial colliding atom pair to an excited pair of 
free atoms and cannot be regarded as photoassociation. 
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IV. COHERENT TRANSIENT DYNAMICS 

Solving the time-dependent Schrodinger equation al- 
lows for a detailed study and analysis of the dynam- 
ics induced by the pump pulse. Figure 0] displays the 
excited state population, J di?|vI'e(_R; t)p, and the fre- 
quency of the induced dipole, i.e. the derivative of the 
time-dependent phase, (/>{p)(t), of J dR'i'*{R; t)fi^e{R\ t)) 
(red dashed lines in Fig. |3]a, c, d, f) as a function of time 
and compares them to the envelope and the instanta- 
neous frequency of the pump pulse (black solid lines in 
Fig.Hla, c, d, f). Also shown is the relative phase between 
the transition dipole and the field, i.e. 4>{^){t) — 4'E(t) 
(solid black line in Fig. |4]b). Characteristic oscillations 
of the excited state population are observed which are 
due to the long tails of the electric field with cut spec- 
trum, cf. Fig , m c. They correspond to the coherent 
transients [29| that have been observed before in the ex- 
citation of atoms with short pulses, see e.g. Refs. [l8l[l9l|. 
Note that in our case, a molecular signal is observed in 
the ionization detection [2]| . 

The transient interaction of the weakly bound molecule 
and the pulse is governed by the induced transition 
dipole, in particular by the relative phase between the 
induced dipole and the field. The dynamics obtained for 
a single spectral cut, i.e. for excitation below the Dl line 
(left-hand side of Fig. [4|), is discussed first. In this case, 
the instantaneous frequency of the laser pulse, shown by 
the black solid line in Fig. [4] a, coincides with the cut- 
off frequency before and after the pulse maximum but 
dips to 12250 cm~^ during the main part of the pulse. 
The depth of the dip is determined by the spectral band- 
width of the pulse. The temporal tails to both sides of 
the main pulse peak lead to significant field amplitude 
that oscillates with the cutoff frequency, cf. black solid 
line in Fig. d] a. During the main peak of the pulse, 
around t = Q, a, dipole is induced by off- resonant excita- 
tion. It is at first driven by the strong electric field and 
follows its oscillation. After the main peak of the pulse 
is over, the coupling is reduced and the dipole oscillates 



FIG. 4: (color online) Details of ultrafast dipole dynamics. 
a-|-d: instantaneous frequency of the pump pulse (solid black 
line), frequency of induced dipole (dashed red line), position 
of the Dl, resp. D2, resonance frequency (dash-dotted line), 
b+e; relative phase between dipole and light field, c-l-f: elec- 
tric field envelope of pump pulse (solid black line), excited 
state population (dashed red line) for the 0j(5s + 5pi/2) state 
(c) and Og (5s -I- 5p3/2) (f). The vertical dotted blue lines in- 
dicate the position where the relative phase crosses (n -f 
n = 0, 1,2. 



with its intrinsic frequency close to the Dl resonance fre- 
quency, cf. red dashed line in Fig. [4] a. This is due to 
the excited state population being dominantly in vibra- 
tional levels close to the dissociation limit, cf. Section 
mil The oscillation of the excited state population goes 
through one period as the relative phase between dipole 
and field changes by 2 tt, see the dotted blue lines and the 
red dashed curve in Fig. |4]c. The dotted blue lines indi- 
cate the times when the relative phase crosses odd integer 
multiples of tt. These times coincide with the minima in 
the oscillation of the excited state population. In the long 
pulse tail, this oscillation results from a beating between 
the transition dipole and the field that is oscillating with 
the cutoff frequency. As the long tail of the pump pulse 
decays (black solid line in Fig.[4]c), the oscillations of the 
transition dipole subside as well. This damps effectively 
the oscillations in the excited state population since the 
energy exchange between pulse and molecule vanishes. 

The dynamics observed for a pulse to which two spec- 
tral cuts are applied, i.e. for excitation below the D2 
resonance frequency (right-hand side of Fig. is over- 
all similar to that obtained for excitation below the Dl 
resonance frequency. However, two spectral cuts lead to 
a distinct beat pattern in the amplitude and phase of 
the electric field which carries over to the dynamical ob- 
servables transition dipole and excited state population. 
The beat pattern is determined by the positions of the 
two cuts relative to the carrier frequency of the field. The 
induced dipole oscillates around the cut position before 
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the main peak of the pulse. This oscillation is of little 
relevance, however, since the amplitude of the induced 
dipole which is related to the excited state population is 
extremely small. It dips down to 12725 cm~^ during the 
main peak of the pulse and oscillates with its intrinsic 
frequency close to the D2 resonance superimposed with 
the fast beat pattern after the main peak is over. While 
the beat patterns are easily resolved in the calculations 
where a time step of 1 fs was employed (the results shown 
in Fig. [?] are not convoluted with a Gaussian represent- 
ing the effect of the probe pulse), they require a probe 
pulse with FWHM of 100 or less and a sufficiently good 
signal-to-noise ratio to be observed experimentally in e.g. 
the excited state population. In the experiment discussed 
here the comparatively long duration of the probe pulse 
does not allow for the observation of this beat pattern. 

The coherent transients are solely due to electronic 
transitions as explained in detail below in Section |Vl 
Their dynamics can therefore also be modelled in the 
framework of dressed states of a two-level system. This 
approach has the advantage of yielding analytic expres- 
sions which describe the coherent interactions causing the 
dipole modulations [3(]| . 



V. COMPARISON OF THEORETICAL AND 
EXPERIMENTAL RESULTS 

A detailed comparison of the quantum dynamical cal- 
culations to the experimental results is presented in the 
following. A typical pump-probe signal recorded in an 
experiment with excitation frequency near the D2 line 
is depicted by the black curves in Fig. [H while the grey 
curves represent the excited state population obtained in 
the calculations after convolution with the probe pulse 
(cf. Sec. |TT]) . For a chosen cutoff position below the D2 
resonance (here 10cm~^), a constant molecular ion sig- 
nal is measured at negative time delays where the probe 
pulse precedes the pump pulse, cf. Fig. [5^. For zero 
delay, the pump and probe pulses coincide in time and 
a peak of about 0.5 ps width is observed. For positive 
time delays where the pump pulse precedes the probe 
pulse, an increased molecular ion signal occurs with char- 
acteristic oscillations that are fully damped after a few 
picoseconds. An excellent agreement between the mea- 
sured pump-probe data and the simulated excited state 
population in the attractive potentials of the 5s -I- 5p3/2 
is observed for positive delays. The finite offset of the 
experimental signal at negative delays is absent in the 
theoretical curve. The calculations are done for a single 
pump pulse, while the experiments operate with a high 
repetition rate, 100 kHz. The signal at negative delays is 
therefore attributed to molecules in the electronic ground 
state that were produced by the previous pump pulse and 
spontaneous emission j2^. This effect is not taken into 
account in the calculations. 

When the atomic resonance is included in the pump 
pulse spectrum. Fig. [SJs, the total pulse energy needs to 



be reduced by 50% with an optical gray filter to avoid 
a too strong perturbation of the ultracold cloud. In this 
higher ion signal compared to that obtained with 
a cut to the red of the resonance frequency is observed 
for positive time delays. The characteristic oscillations at 
positive delay disappear almost completely, cf. Fig. \Ejp. 
When a transform-limited pulse is considered in the cal- 
culation, the population is quickly transferred to the ex- 
cited state, and no oscillations are observed in Fig. [5}: (the 
inset shows a zoom of 200 fs around zero delay). In the 
experiment, a transform- limited pulse led to a destruc- 
tion of the sample. Fig. [5] illustrates that the transient 
oscillations depend strongly on the spectral cutoff of the 
pump pulse spectrum. 

Fig. [5] shows the time- resolved molecular ion signal for 
different spectral cutoffs below the Dl and D2 line res- 
onances in comparison with the corresponding quantum 
dynamical calculations. Generally, the pump-probe spec- 
tra obtained by exciting the ultracold ensemble near the 
Dl and D2 line resonances are very similar. A notable 
difference is observed in the RbJ count rate which is 
larger for excitation near the D2 line and positive delays. 
This is attributed to the larger dipole moment for the 
5s — > 5p3/2 transition. When the sample is irradiated by 
the probe pulse only, a constant level of molecules is de- 
tected. This signal is due to molecules that are produced 
by the trapping light. It is subtracted as background 
from each pump-probe spectrum. 

In Fig. [S] the 'detuning' of the spectral cut position 
with respect to the atomic resonance increases from top 
to bottom. A decreasing mean ion signal at positive time 
delays is observed for increasing cutoff position for exci- 
tation close to both the Dl and D2 line resonance. More- 
over, the characteristic oscillations at positive time delay 
decrease in their period and amplitude as the cutoff po- 
sition is moved further away from the resonance. The 
decrease in oscillation period is explained in terms of the 
instantaneous frequency of the pump pulse. A linear de- 
pendence of the oscillation frequency on the frequencyof 
the cutoff position was recorded, cf. Fig. 2 of Ref. [20| . 
This is due to the laser field oscillating with an instanta- 
neous frequency that corresponds to the cutoff position, 
cf. Sec. IIVI The decrease in amplitude with increasing 
cutoff detuning is attributed to the most important fre- 
quencies being excluded from the pulse. Moreover, the 
pulse energy is reduced when more frequency components 
are removed from the pump pulse spectrum. Then the 
instantaneous intensity and hence the coupling between 
the field and and the molecular dipole become smaller. 

One might expect nuclear dynamics of the excited state 
wavepacket to show up as oscillations in the molecular 
ion signal. However, the vibrational periods of the ex- 
cited state levels that are populated, cf. Fig. O are at 
least one order of magnitude larger than the oscillation 
periods observed in Fig. [6l Moreover, one would expect 
to see wave packet revivals in the vibrational dynamics 
which have not been observed even when scanning the 
delay for very long times (up to 250 ps). In the calcu- 
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FIG. 5: Experimental (black) and theoretical (grey) pump-probe signals for pump pulses with spectral cutoff (a) 10cm~^ red 
to the D2 line resonance, (b) directly on the resonance frequency, and (c) theoretical signal for a pump pulse without spectral 
cutoff (transform- limited pulse). 



lations, vibrational dynamics should yield slightly differ- 
ent oscillation periods for the various potentials due to 
the different effective C3 coefficients. However, identical 
oscillation periods were obtained for all potentials, inde- 
pendent of their character. All these observations clearly 
rule out molecular vibrational dynamics as origin of the 
observed oscillations and support their purely electronic 
character. 

Fig. [7] displays the variation of the averaged asymptotic 
pump-probe signals with the cutoff position. Asymptotic 
refers to delay times for positive delay {t > 0) where the 
oscillatory behavior is fully damped and to delay times 
before the signal rise for negative delay {t < 0). The 
measured data are shown as upward triangles for pos- 
itive delay and as downward triangles for negative de- 
lay. The results of the calculation were scaled to match 
the experimental data and are represented by a grey line 
(positive delay only) . The insets display the ratio of the 
asymptotic pump-probe signals at positive and negative 
delay versus the cutoff position. This ratio is obtained 
to be about 0.3 at the Dl line and about 0.1 at the D2 
line. As discussed in detail in Ref. [2l|, the nearly con- 
stant ratio indicates that the signal at negative delay, 
when the probe precedes the pump pulse is related to 
molecules that were photoassociated by a previous pump 
pulse and have undergone spontaneous decay to the elec- 
tronic ground state. The experimental result for positive 
delay is well reproduced by the calculation. In particu- 
lar, in accordance with the results of Fig. [S] an increase 
of the excited state population is observed as the cutoff 
position approaches the resonance frequency. 

The influence of a linear frequency chirp on the pump- 
probe signals is discussed with the help of Fig. [S] 
Measured (calculated) pump-probe spectra are shown for 
different positive chirps in Fig. [8^ (c) and for negative 
chirps in Fig. [5)3 (d) for an excitation near the D2 line 
resonance and a cutoff position of 8cm~^ below the res- 
onance frequency. The spectral cutoff in combination 
with the applied quadratic phase shift causes strong dis- 
tortions of amplitude and frequency behavior during the 
pulse compared to the transform-limited and spectrally 
cut, but unchirped pulses. A frequency chirp conserves 



the spectral bandwidth but prolongs the pulse. Since the 
pulse energy remains unchanged, the peak intensities are 
decreased. This leads to a reduction of the detected ion 
signal compared to an unchirped pump pulse for both 
negative and positive time delays, cf. the black and grey 
lines in Fig.[8l Coherent transients are present also in the 
pump-probe spectra for chirped pump pulses. The chirp 
does not change their period, but it induces a shift of 
the zero delay peak, as indicated by the arrows in Fig. [S] 
This behavior is well reproduced by quantum dynamical 
calculations. For a positive (negative) chirp, the larger 
instantaneous frequencies occur after (before) the smaller 
ones during the pulse. The larger instantaneous frequen- 
cies are more likely to cause the off-resonant excitation of 
population into the excited state than the smaller ones. 
This explains the shift of the zero delay peak toward later 
(earlier) times for positive (negative) chirps. 

In previous calculations with chirped pulses, e.g. 0, Q 
an enhanced efficiency was proposed for chirped pulses. 
Here, an enhancement of the excited state population 
by using linearly chirped pump pulses was not observed. 
This is due to the off-resonant character of the excita- 
tion, cf. Sec. mil References ^ ^] describe resonant 
excitation, i.e. only levels with free-bound transition fre- 
quencies within the spectrum of the pulse were addressed. 
Moreover, picosecond pulses were considered where some 
vibrational dynamics occurs during the pulse. The pro- 
posed enhancement then results from an interplay be- 
tween the wavepacket rolling down the potential and 
the coupling with the field. This cannot be expected 
in our case. Note that such an interplay between nuclear 
and electronic dynamics has been recently observed for 
nanosecond pulses, demonstrating coherent control of ul- 
tracold collisions with frequency-chirped light [31] ■ In ad- 
dition to reducing excitation efficiency, the linear chirps 
affect the dynamics of the induced dipole. This can be 
seen in the shift of the first peak relative to the oscilla- 
tions indicated by the arrows in Fig. [5] for positive and 
negative chirps. As discussed above, the oscillations arise 
from the energy exchange between the induced dipole and 
the electric field of the pulse. They can therefore be ex- 
pected to be very sensitive to phase manipulations of the 
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FIG. 6: Experimental (black) and theoretical (gray) pump-probe signals for pump pulses with different spectral cutoff positions 
red to the Df line resonance (5s + 5pi/2-asymptote, left) and D2 line resonance (5s + 5p3/2-a'Symptote, right), respectively. The 
constant level of molecules measured by irradiating with the probe pulse only is indicated by the dashed lines, and the zero of 
the theoretical results is shifted to coincide with this level. 



pulse which also affect the relative phase between dipole 
and field. 

Since the molecular dynamics is reduced to the elec- 
tronic dynamics, i.e. to that of a two-level system, the 
interaction between the molecule and the field can also 
be analyzed in terms of the evolution of the Bloch vector 
[soi l . Fig. [9] shows the dispersive and absorptive compo- 
nents, U{t) = Pge{t)+Pegit) and V{t) = i{pge{t) - Peg{t)), 

of the Bloch vectors for chirped and unchirped pulses ver- 
sus time (s^ j. For unchirped pulses, cf. middle graph in 
Fig. [HI the Bloch vector moves merely through a single 
loop in the u — v plane during the main pulse peak. For 
chirped pulses the Bloch vectors perform rapid rotations 



in the u — v plane during the main pulse peak as the 
instantaneous frequency of the pulse varies, cf. Fig. 5^. 
For negative (positive) chirps, cf. the left (right) graph 
of Fig. [HI the rapid rotations occur before (after) t = 0. 
As a result the final position of the Bloch vector in the 
u — v plane is different for different sign of the chirp. The 
modulations on the absorptive component that are re- 
sponsible for the energy exchange between molecule and 
field then oscillate with a phase offset compared to zero 
chirp. 
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the ratio of the asymptotic signals at positive and negative delay for a given cutoff position. 



VI. CONCLUSIONS 

We have presented theoretical calculations describing 
femtosecond pump-probe photoassociation experiments 
[20l . [2]| . Two colliding rubidium atoms in their elec- 
tronic ground state are excited by the pump pulse into 
bound levels and continuum states near the 5s -I- 5pi/2 
and 5s -|- 5p3/2 asymptotes. They are ionized by a sec- 
ond, time-delayed probe pulse and detected by a mass 
spectrometer. The model considers only the pump pulse 
explicitly while accounting for the finite width of the 
probe pulse by convolution with a Gaussian. The pump 
pulse was modeled according to the experimental con- 
ditions, fn particular, spectral cuts of amplitude and 
quadratic spectral phases were considered. Photoasso- 
ciation is observed as excitation of the two atoms into 
weakly bound molecular levels, mostly very close to the 
dissociation limit. It is attributed mainly to off-resonant 
excitation. The spectral cutoff produces pulses in time 
with very long tails giving rise to transient oscillations 
of the ion signal, respectively excited state population. 
Coherent transients were analyzed in terms of the field- 
induced dipole and the relative phase between system 
and field. Pump pulses obtained for different cut posi- 
tions and linear chirps were employed. The oscillation 
period of the coherent transients is directly related to 
the cutoff position. A negative (positive) chirp simply 
shifts the zero peak delay to shorter (larger) delay times. 
Calculations were performed for a number of different 
excited state potential energy curves correlating to both 
the 5s -I- 5pi/2 and 5s -I- 5p3/2 dissociation limits yield- 
ing, however, identical oscillation periods of the coherent 
transients. This indicates a purely electronic character 
of the observed oscillations. The analysis of the two-level 
dynamics was further confirmed by inspection of the sys- 
tem's Bloch vector. 

The experiments modelled in this paper have provided 



the first evidence for femtosecond photoassociation. The 
off-resonant excitation of molecular levels and the coher- 
ent transients both clearly demonstrate that the observed 
phenomena are beyond the sudden approximation. How- 
ever, the molecules created by the femtosecond pump 
pulse are extremely weakly bound - they do not provide a 
suitable starting point for the application of a dump pulse 
in order to create ultracold molecules in their electronic 
ground state. In future experiments, the pulses should 
therefore be better adapted to the timescales of the prob- 
lem. In particular, a bandwidth of a few wavenumbers 
should be best to create a spatially localized excited state 
wavepacket [si] that can subsequently be dumped to the 
electronic ground state [6, 7J. 



Theoretical estimates predict a moderate photoasso- 
ciation efficiency of about one molecule per pulse for a 
magneto-optical trap (MOT) [8] . The efficiency is limited 
by the comparatively small pair density at the distances 
where photoassociation is effective. Future experiments 
on short-pulse photoassociation should therefore be com- 
bined with efforts to increase the pair density. Two routes 
can be pursued towards this goal, (i) The atoms can be 
cooled further to reach the quantum degenerate regime. 
Compared to the highest possible densities in a MOT, 
about three orders of magnitude can be gained this way. 
(ii) External fields can be employed to enhance the pair 
density at short and intermediate distances. For ex- 
ample, an increase in the photoassociation efficiency of 
up to five orders of magnitude is predicted for magnetic 
field control in Feshbach-optimized photoassociation fs^ . 
Pursuing one of these routes or a combination thereof, an 
efficiency of 10"^ molecules per pulse or more is to be ex- 
pected in photoassociation with picosecond pulses. 
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FIG. 8: (color online) Pump-probe spectra for different quadratic spectral phase shifts, i.e. linear chirps and a pump pulse 
cutoff position of 8cm~^ below the D2 line resonance frequency. Experimental (theoretical) results are shown in top (bottom) 
panels. As indicated by the arrows, a positive chirp (a, c) induces a shift of the zero delay peak to positive time delays, whereas 
a negative chirp (b, d) causes the opposite effect. 
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